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ABSTRACT 


We report on monitoring observations of the TeV 7 -ray binary 
HESS J0632-I-057, which were carried out to constrain the interaction between 
the Be circumstellar disk and the compact object of unknown nature, and provide 
for the hrst time high-dispersion (R> 50 000) optical spectra in the second half 
of the orbital cycle, from apastron through periastron. The Ha, H/9, and H 7 line 
prohles are found to exhibit remarkable short-term variability for ~1 month after 
the apastron (phase 0 . 6 - 0 .7), whereas they show little variation near the perias¬ 
tron. These emission lines show “S-shaped” variations with timescale of ~150 
days, which is about twice that reported previously. In contrast to the Balmer 
lines, no prohle variability is seen in any Fe II emission line. We estimate the 
radii of emitting regions of the Ha, H/9, H 7 , and Fe II emission lines to be ~30, 
11, 7, and 2 stellar radii (R*), respectively. The amplitudes of the line prohle 
variations in different lines indicate that the interaction with the compact object 
affects the Be disk down to, at least, the radius of 7 R* after the apastron. This 
fact, together with little prohle variability near the periastron, rules out the tidal 
force as the major cause of disk variability. Although this leaves the pulsar wind 
as the most likely candidate mechanism for disk variations, understanding the 
details of the interaction, particularly the mechanism for causing a large-scale 
disk disturbance after the apastron, remains an open question. 

Subject headings: binaries (including multiple); close — gamma rays: stars — stars: 
emission-line. Be — stars: individual: (HESS J0632-I-057) — X-rays: binaries 
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Introduction 


TeV 7 -ray binaries are a subclass of binaries with a compact object, established in 
the 2000’s (e.g., Dubus 2013 


, for a recent review). These systems have a spectral energy 


distribution with a peak beyond 1 MeV and are variable at multi-wavelengths up to TeV 
energies. There are 5 known binaries of this kind, all of which have either an 0-type 
main-sequence star or a Be star with a circumstellar disk as the optical counterpart. The 
nature of the compact object is unknown in all systems but one (PSR B1259-63). For 
such systems two competing scenarios have been proposed, based on the different nature 
of the compact object, and hence the different mechanism of high energy emissions. One 
scenario assumes that the collision between a relativistic pulsar wind and a stellar wind 


and/ or circumstellar dis 


(e.g.. 


c prod uces strong shocks, where the high energy emission arises 


Maraschi fc Treves 


198l[) . This pulsar wind scenario has successfully been applied to 


PSR B1259-63. On the other hand, the other scenario assumes the presence of an accreting 
black hole (or neutron star). In this microquasar or accretion/ejection scenario, a large 
amount of mass transfer red from a companion star p owers relativistic jets, where y-ray 


emission originates (e.g.. 


Levinson fc Blandford 


19961). 


HE SS J0632+057 f06^32" ^59h25 + 05°48'01".2) is a recently established TeV y-ray 


binary (lAharonian et ah 


20071) comprised of a BOVe stm and a compact object of unknown 


nature. The orbit is wide (Porb =31514 days; 


Alin et ah 


Casares et al. 


20121). T 


bands flAliu et al. 


2014J ) and highly eccentric (e =0.83; 


le system shows correlated variability in the X-ray and TeV energy 


2014f) . with two peaks in one orbital cycle: the primary outburst prior 


to apastron (orbital phase 0 = 0.3-0.4) and the secondary outburst after it (0 ~ 0.6-0.9), 
with an X-ray dip in between. This is puzzling, given that these outbursts and the dip 
occur when the compact object is orbiting very far from the Be star. Understanding the 
mechanism for these phenomena is one of the key issues of this system. 
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The optical counterpart of HESS J0632+057, MWC 148, is established to be a classical 
Be star, i.e., a massive star wi th a geometrically th i n, circumstellar disk, w here Balmer and 


other emission lines originate flRivinius et ah 


20181 


Aragona et ah! (120101) observed the Ha 


emission line after apastron (0 ~ 0.7-0.8). Their p rohles over continuou s 35 days showed an 


‘S-shaped” variability with a period of ~ 60 days. 


Casares et af 


(120121 ) reported the orbital 


modulation in the Ha line prohle parameters, likely caused by the interaction between the 
Be disk and the compact object. However, the lack of observations in 0=0.6-0.8 makes it 
difficult to constrain the modulation after apastron. 


In order to constrain the interaction between the Be disk and the compact object, 
we have been monitoring HESS J0632+057 using various methods such as high-dispersion 
spectroscopy, photometry, and polarimetry. In this paper we report the initial result of 
the line prohle variabilities over 160 days from 2013 October to 2014 April, covering the 
second half of the orbital period. In Sect. 2, we summarize our observations. The results 
are described in Sect. 3. In Sect. 4 we discuss the Be disk region affected by the compact 
object and the nature of the interaction. 


Observations 


Optical high-dispersion spectroscopic observations of HESS J0632-I-057 were carried 
out on 12 nights from 2013 October 31 (phase 0.555) to 2014 March 11 (phase 0.969) at the 
Okayam a Astrophysical Ob servatory with a 188 cm telescope and HIDES with hber-feed 


system (IKambe et ah 


period of 315 days from 


20131). Here , we ca lculated the orbital phase, 0, takin g the orbital 


Alin et ah 


( 20141 ) and the or igin (JD 2454857.5 


(120121) ■ and setting the periastron phase to be 0 (In 


Note that, although 


Casares et ah 


2012 


from 


Casares et ah 


Casares et ah 


, it is set to 0.967). 


(120121) determined the ephemeris usi ng the orbi 


of 321 days, the orbital parameters remain approximately unchanged (lAliu et al 


al pe riod 


2014h . 
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The observed spectra covered 4200-7400 A wavelength range, with S/N of ~ 150 aronnd 
Ha. The typical wavelength resolution i? is ~ 50 000. The data were reduced in the 
standard way, using the IRAlji] echelle package. 

High-dispersion spectra were also taken on two nights (2014 February 05, phase 
0.861, and April 10, pha se 0.064) using the Canada France Hawaii Telescope/ESPaDOnS 
( Manset fc Donati LoOS ) in spectropolarimetric mode. The spectra covered a wavelength 
range of 3700-10 500 A with a resolving power of ~ 68 000. We obtained reduced data 
using the Libre-Esprit/Upenay pipeline, provided by the instrument team. We rectified the 
normalized intensity by re-determining the continuum level around each line, in order to 
compare with HIDES data. In this paper, we focus on the spectroscopic variability. The 
polarimetric data will be discussed in a forthcoming paper. 


In addition to spectroscopic data, V-, Re-, and Jc-banc 


obtained from 2013 September to 2014 May using HOWPol flKawabata et al 


HONIR f Sakimoto et al. 


photometric data we re 


2008) and 


20121 ) attached to Kanata 1.5m Telescope at Higashi-Hiroshima 


Observatory in Hiroshima, Japan. We observed HESS J0632-I-057 on 73, 71 and 10 nights 
in V, Rc, and Ic bands, respectively, using HOWPol. We used HONIR to observe the 
source on 8 nights, in these three filters. The logs of the photometric and spectroscopic 
observations are given in Table [H 


3. Results 

Figure [T] shows the averaged Ha, H/?, and Hy profiles. In these profiles, the vertical 
scale of the Ha line prohle is half that of the other line prohles, because the flux is 

^http://ir af.noao.edu/ 

^http://www.cflit.hawaii.edu/Instruments/Upena/ 
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much stronger in the Ha. The Ha line is very strong [EW ~ —30 A) and signihcantly 
asymmetric. The prohle exhibits a double peak at ~ 0 km s“^ and ~ 70 km s“^, with the 
latter being brighter than the former, and a hump at ~ —100 km s“^ on the blue wing. On 
the other hand, the averaged H/3 line prohle {EW ~ —3.9 A) is rather symmetric, with a 
slightly stronger blue peak. The Hy line, which is on a broad absorption component, also 
exhibits a double-peaked prohle with stronger blue peak. 


Balmer lines exhibited complicated line prohle variabilities during the monitoring 
period. Figure [2] displays all observed prohles of the Ha, H/?, and Hy lines, while Fig. E] 
presents the time sequence of the residual spectra of these lines from the average. The 
variation in the eqnivalent width, EW, the fnll width at half maximnm, FWHM, and the 
centroid velocity, 14, of the Ha line prohle are shown in Fig. 01 Note that FWHM and 14 
are measured by htting the whole prohle with a single Gaussian and the line wings with a 
Voigt prohle, respectively. 


In the long term (> 100 days, from 0 ~ 0.55 — 1), the Ha line changed from a 
red-enhanced prohle to a rather symmetric prohle, while the stronger peak of the H/? and 
Hy lines changed from the red side to the bine side. In the residnal spectra of these lines, 
S-shaped variability is visible; the bright part at ~ 200 km s“^ migrates from the bine side 
to the red side at 0 ~ 0.7, and returns to the blue side at 0 ~ 1.0. Simultaneously, the 
faint part at ~ —200 km s“^ migrates in the opposite direction. At the same time, the 
maximnm and minimnm, respectively at ~ 100 km s“^ and ~ —100 km s“^ in the residual 
spectra, migrate in a similar fashion. The period of these variations is estimated at ~ 150 
days using Fonrier an alysis, which is abont twice the ~ 60-day variation period reported by 


Araeona et ah 


(120101 ). We will discuss this in more detail in the next section. On the other 
hand, in the short term (<C 100 days), a bright hump appears at ~ 100 km s“^ in Balmer 


line prohles at 0 ~ 0.6. This hnmp disappears by 0 = 0.65 (see the third, forth, and hfth 
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profiles from the bottom in Figs. [2]and|3]). When the hump appears, EW(Ha) increases by 
~ lA and FWHM(Hq!) decreases by ~ 10 km s“^. These long- and short-term variabilities 
respectively seem to occur in phase in different lines, although time lags between lines 
cannot be ruled out because of the low cadence of the observations. 


Recently, 


Casares et ah 


(I 2 OI 2 I) reported the orbital modulation of the profile parameters 


of the Ha line. The ir data cover the phas e intervals of 0-0.6 and 0.8-1, using the orbital 


period of 321 days flBongiorno et ah 


201 ll) . If we t ake the same o r bital 


spans the orbital phase 0.4-1, which overlaps with 


Casares et ah 


period, our data 


I 2 OI 2 I ) in the phase 


range 0.4-0.6 and 0.8-1 and provides for the hrst time the modulation data in the phase 
interval of 0.6-0.8. In Fig. HJ EW(Ha) seems to increase (phase 0.4-0.5) an d then decrease 


back 


;o the previous value in about 30 days (phase 0.5-0.6), as reported by 


Casares et ah 


(120121 ). Afterwards, it seems to fluctuate for about 60 days (phase 0.6-0.8), which could 


be associated with the variabilities mentioned above. Vc{lia), whose var iation is thoug ht 


to refl ect the orbital motion of the Be star, shows a similar pattern to 


Casares et af 


(120121 ) except just after the apastron and just before the periastron. This difference in 


Vc is possibly caused by the profil e variabilities. FWH MfHaj, on the other hand, shows 


a different variation pattern from 


Casares et ah 


(120121) . It stays constant ex cept around 


apastr on (phase 0.45-0.55 for the period of 321 days) in this work, whereas in 
(I 2 OI 2 I) it showed a sinusoidal pattern. 


Casares et ah 


Although many Fe II lines are contaminated by neighboring lines, there are six 
less-contaminated lines (A =5018, 5316, 5363, 5535, 6433, and 6456 A) that enable us to 
analyze the line profile variabilities. In our analysis, we have used these six lines. The 
averaged Fe II A5363 line profile in Fig. [1] is shown as a representative Fe II profile. 

In contrast to the Balmer lines, the Fe II emission lines show symmetric double-peaked 
prohles, as predicted for a rotating, axi-symmetric disk. Figures [2] and [3] also display the 




























time sequence of the observed and residual Fe II A5363 pro f iles, r espectively. As seen in 


these hgures, Fe II lines exhibited no variation. 


Arias et al 


(120061) analyzed the spectra of 


several Be stars and conhrmed that Fe II emission lines arise from a disk region of radius 
of 2.0 ± 0.8 R*. They also found empirical relationships between the projected rotational 
velocity of the central star [V sinz) and the prohle parameters. Applying their relationships, 
we have derived V sinz to be ~ 230 — 240 km s“^. 


Finally, we estimate the radii of emitting regions of Balmer lines. For the Ha line, 
we use the mean peak separation of the last six prohles (0 > 0.8), because the averaged 
prohle in Fig. [H in which all our observed epochs were used, has too complicated features 
to determine the peak velocities. For the H/3 and Hy lines, we use the averaged prohles. 
As a result, we have the peak separations of the Ha, H/?, and Hy lines to be 90 km s“^, 
143 km s“^ and 176 km s“^, respectively. Adopting V sinz derived above, we obtain the 
radii of the emitting region of these lines as ~ 30 R* (Ha), HR* (H/5), and 7 R* (Hy), 
where R* is the radius of the Be star. 

Throughout our monitoring period (from apastron to periastron), optical brightness of 
HESS J0632+057 stayed constant in the range of 9.05-9.32 {V) mag, 8.53-8.71 (Rc) mag, 
and 8.58-8.67 (R) mag, within the typical 1-a error of 0.05 mag. The average values are 
V = 9.13 ± 0.05 mag, Rc = 8.58 ± 0.06 mag and Ic = 8.62 ± 0.03 mag. 


4. Discussion 

4.1. The Disk Radius Affected by the Compact Object 

Although the monitoring was performed mostly after apastron, the line prohles 
exhibited remarkable variabilities. Wide wavelength coverage revealed that the variations 
are seen not only in the Ha line but also in the H/5 and Hy lines. No variation in Fe II lines 
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indicates the inner part of the Be disk was kept undisturbed during the observing period. 


This fact is in agreement with little variability in the optical brf 


to originate from the disk region within 2-3 R* flRivinius et ah 


ghtness, which is thought 


20131). The slight, but 


significant, variation in the Hy line profile implies that the interaction with the compact 
object affects the Be disk down to, at least, a radius of 7 R*. 


4.2. S-shaped Variation 


T 

E]), as 


le dynamical r e sidua l spectra of the Balmer lines showed an S-shaped variation (Fig. 


Araeona et ah 


(I 2 OIOI) reported. High-dispersion residual spectra clearly show that 
there are two pairs of peaks at ~ —200/200 and ^ —100/100 km s“^. The lower velocity 


peaks are the same as those of 


Aragona et al 


1 2010 !) (see their Fig. 3). The 


varia tion is ~150 days, which is about twice the ~ 60-day variation period flAragona et al 


period of the 


2OI0I). 


The S-shaped variation seen in many Be stars is thought to be caused by global disk 


oscillations. T 


ley are mostly low-frequency tu = 1 o scillations, where m is the azimuthal 


wave number fjOkazaki 


1991 


Papaloizon et al. 


19921 ). but in eccentric binaries, it is also 


possible that m > 1 osc illations are excited by the 


the tidal potential (e.g.. 


Artvmowicz fc Lubow 


corresponding Fourier components of 


19941) . If the period of observed variation is 


~ 150 days, about half the orbital period, it might be due to an m = 2 oscillation mode in 
the Be disk, which could be excited by the m = 2 Fourier component of the tidal potential 
of the compact object. 


The difference of the periods between this work and 


Aragona et al. 


( 2010 ) might be 


caused by the difference of the length and frequency of the observations. Because we 
observed HESS J0632-I-057 for about half an orbital cycle (~ 160 days), variabilities of this 
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timescale can be detected. With low cadence of observation, typically once every 10-20 
days, however, prohle vari abilities of shorter tim escales are invisible. On the contrary, the 


short monitoring epoch of 


Aragona et ah 


2OIOI) (35 days) wonld have hidden the variability 


with 


onger timescales (> 100 days). Moreover, observations of this work and 


Aragona et ah 


(12010) have an interval of ~1800 days. Therefore, there are possibilities that the oscillation 


mode has changed or that both oscillation modes exist. 


4.3. Short-Term Variability 


In addition to the S-shaped variation, a remarkable variability was seen in the Balmer 
lines after apastron ((;/)=0.60-0.65). The short lifetime of this variability (< 50 days) 
indicates that it is temporarily cansed by an external force snch as the tidal force and 


the ram pressnre of the 
lifetimes ( 1000 days; 


pulsar wind. 


Okazaki 


1991 


Decause internal waves/oscillations have much longer 
and references therein). It is surprising, however. 


that such a remarkable variability appeared after apastron, when the interaction is expected 
to be weak because of the large distance between the Be disk and the compact object. In 
order to understand this phenomenon, it is essential to constrain the starting phase and 
repeatability of the variability. 


Comparison of the EW(Ha) between this work and 


Casares et af 


( 2012 ) implies the 


presence of a regular orbital modulation, but different absolute values due to different 
spectral resolutions make it difficult to quantitatively analyze it. Given the orbital period 
close to one year, it is necessary to monitor over two or three more successive cycles, within 
the timescale of Be disk variability (< 1000 days), in order to cover the full orbital phase 
with the same spectral resolution. 
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4.4. Implications for the Nature of the System 

Because the compact object is located far away from the Be star around apastron 
(~ 100 R*), the gas in the Be disk at the radius of 7 R* cannot be signihcantly affected by 
the tidal force of the compact object, whose strength is less than 10“^ of the gravity of the 
Be central star at this radius. This leaves the interaction with the pulsar wind as the most 
likely mechanism for the Balmer line variability. 

At the periastron passage, on the other hand, the compact object is thought to pass 
through the Be disk at ~ 10 R*. Nonetheless, the H/5 line, emitted from the disk region of 
the similar radius, showed no remarkable variation at 0 =0.06 (see the top profile in the 
upper right panel of Fig. H]). This fact suggests that either the tidal interaction is very weak 
in this system or the orbital period is longer than 315 days, so that the compact object had 
not passed the disk yet at the time of observation. The former possibility could be realized 
if the Be disk rotates in the retrograde direction or is tilted by a large angle with respect 
to the orbital plane, where the disk gas and the compact object interact on a very short 
timescale. Future observations at significantly later phases will distinguish between these 
two possibilities. 


4.5. Flip-Flop Scenario 


As described in Sects. 14.31 and 14.41 our observations leave the pulsar wind model as 
the sole candidate for HESS J0632+057. If the compact object is a pulsar, little line profile 
variation around periastron (the top two profiles in each panel of Fig. |3]) sugg ests that the 


pulsar wind is too weak to significantly affect the Be disk during these phases. 


Torres et ah 


(120121 ) proposed a fiip-fiop scenario for another gamma-ray binary, LS I-|-61°303, where the 
pulsar is in a rotationally powered regime in the apastron, while it is in a propeller regime 
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in the periastron (see also iPapitto et al.l l2012l) . In a flip-flop system, if the gas pressure of 
the Be disk overcomes the pulsar-wind ram pressure, the pulsar wind is quenched. Because 
the Be disk of HESS J0632-I-057 is estimated to be about three times larger than the binary 
separation at periastron, the compact object crosses a dense region of the disk near the 
periastron. In such a situation, the strong gas pressure is likely to quench the pulsar wind 
and hence suppress high-energy emissions. 

In the framework of the pulsar wind model, there are a few mechanisms that might 
explain the short-term episodic variability discussed above. The Be disk is likely as 
l arge as the binary orbit bec ause of no tidal truncation in highly eccentric, large orbit 


flOkazaki fc Neguerue 


Carciofl fc Biorkman 


a 


20011) . Because the disk density rapidly decreases with radius (e.g.. 


20061) . the wind from the pulsar close to the outer part of the disk 


effectively changes its structure, giving rise to remarkable variability in emission lines arising 
from the disk outer part. If the Be disk is misaligned with the orbital plane and a node 
happens to be in the direction corresponding to the phase of the variation (0 ~ 0.6 — 0.65), 
the pulsar-wind effect on the Be disk will be the strongest in this phase interval. After 
0 ~ 0.8 the pulsar comes close to the denser disk region, which is not easily affected by the 
pulsar wind. This might be a cause of the short-term episodic variation after aspastron. It 
is not clear, however, how the pulsar wind can affect the inner disk at the radius of 7 R*. 
Alternatively, the short-term post-apastron variations might be explained as the emission 
from the gas captured by the pulsar, if the pulsar wind is not strong enough to expel the 
surrounding gas. This picture seems to fit well with the flip-flop model. At any rate, 
observational investigation in the first half of the orbital cycle is needed in order to further 
test the pulsar-wind scenario. 

The next periastron pass age of HES S J063 2-I-057 will take place in 2015 December, 


according to the ephemeris of 


Alin et ah 


(120141) (2016 January if the ephemeris is taken 
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from 


Casares et al. 


20121) . Observations covering this period will provide more clnes to the 


complex interaction and the natnre of the compact object in this pnzzling 7 -ray binary. 


This paper is based on the observations taken at the Okayama Astrophysical 
Observatory. We are gratefnl for Dr. Eiji Kambe to kindly observe HESS J0632-I-057. 
This work was snpported by Research Fellowships for the Promotion of Science for Yonng 
Scientists (YM, KT). ATO acknowledges snpport by the JSPS Grant-in-Aid for Scientific 
Research (24540235) and a research grant from Hokkai-Gaknen Edncational Fonndation. 
AGG acknowledges snpport from GNPq (grant 307076/2012-1). 
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Table 1: Observation log. 


spectroscopic data 


Date 

JD 

phase 

OAO/HIDES 





2013.10.31 

2456597.278 

0.555 


2013.11.08 

2456605.275 

0.580 


2013.11.11 

2456608.274 

0.589 


2013.11.19 

2456616.199 

0.615 


2013.12.01 

2456628.219 

0.653 


2013.12.19 

2456646.201 

0.710 


2013.12.29 

2456656.007 

0.741 


2014.01.02 

2456660.047 

0.754 


2014.01.15 

2456673.095 

0.795 


2014.01.23 

2456681.094 

0.821 


2014.02.08 

2456696.968 

0.871 


2014.03.11 

2456728.008 

0.969 

CFHT/ESPaDOnS 





2014.02.05 

2456693.843 

0.861 


2014.04.10 

2456757.772 

0.064 

photometric data 


Filter 

JD 

nights 

Kanata/HOWPol 





V 

2456551-2456795 

73 


Rc 

2456551-2456795 

71 


Ic 

2456722-2456795 

10 

Kanata/HONIR 





V, Rc, Ic 

2456743-2456769 

8 
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Fig. 1.— Averaged profiles of Ha, H/9, H 7 , and Fe II A5363 in the velocity reference frame. 
For illustration purpose, there are vertical offsets between each prohle, and the Ha line flux 
3 V 0.5 . The vertical dashed line indicates the systemic velocity (48.3 kms“^; 


is multiplied 

3V 0.5 

Casares et ah 

2012 ) 
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Fig. 2.— Observed Ha, H/3, H 7 , and Fell A5363 line profiles. Each profile is vertically 
shifted according to the orbital phase, labeled on the right axis. The vertical line is the same 

as Fig. [H 
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Fig. 3.— Dynamical residual spectra of Ha, H/5, Hy, and Fe II A5363 from the average 
profile. The vertical line is the same as Fig. [1] 
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Fig. 4.— Variations of EW (top), FWHM (middle) and Vc (bottom) of the Ha line profile. 
The orbital phase is annotated on the top of the figure, for two different orbital periods 
(Porb) of 315 and 321 days. The vertical solid and dashed lines mark the periastron for the 
former and latter peri ods, respectively. T he horizontal line in the bottom panel indicates 


the systemic velocity flCasares et ah 


2OI2I). 












